The spin resonance and high frequency optical properties of the cuprates. 
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We argue that recently observed superconductivity-induced blue shift of the plasma frequency 
5u)pi in Bi2Sr2CaCu20s+s CI is related to the change in the integrated dynamical structure factor 
associated with the development of the spin resonance below Tc- We show that the magnitude of 
Suipi is consistent with the small integrated spectral weight of the resonance, and its temperature 
dependences closely follow that of the spin resonance peak. 

PACS numbers: 71.10.Ca,74.20.Fg,74.25.-q 



The importance of the resonance spin mode for the 
physics of the cuprates continue to be the subject of in- 
tensive debate. In a generic superconductor, the pair- 
ing of fermions drastically reduces the damping of collec- 
tive spin degrees of freedom at energies below 2 A. For 
a d— wave superconductor, the residual interaction be- 
tween spin fluctuations and fermions gives rise to the ad- 
ditional effect - the development of the exciton mode 
below 2 A (see e.g. Ref. ^ and references therein). 
This mode exists for bosonic momenta near (tt, tt) and 
is commonly called the "spin resonance" . It has been ob- 
served in three different families of high- Tc superconduc- 
tors 0iSISIEI3- This mode is not a "glue" to supercon- 
ductivity as it emerges only in the superconducting state 
(more precisely, below the pseudogap temperature), but 
it affects electronic properties of the cuprates in the su- 
perconducting state. 

Much of recent works on the effect of the spin reso- 
nance on electrons was concentrated on whether the in- 
teraction with the resonance is capable to explain exper- 
imentally detected low- energy features in the fermionic 
spectral functioii, tunneling density of states and optical 
conductivity 0, Ig ■ An example of this behavior is the 
peak/dip/hump structure of the spectral function j^. 

The spin-resonance scenario for the cuprates was re- 
cently questioned lOj on the basis that the measured 
spectral weigh of the resonance is only 1 — 2 percent of the 
total magnetic _spectral weight. Ar. Abanov et al., how- 
ever, argued jTj that this smallness merely reflects the 
fact that the resonance exists only in a narrow momen- 
tum range near Q = (7r/a,7r/a) [roughly, between Q and 
0.8Q], while at momenta where the resonance does exist, 
it strongly couples to fermions. They further argued that 
at small uj ~ 2A, the fermionic self-energy change be- 
tween normal and superconducting states chiefly comes 
from spin fluctuations with energies comparable to w, 
and from spin momenta |q— Q| co/vf, where Wi? is the 
bare Fermi velocity. As vp ie.V * a is large, typical 
spin momenta for uj < lOOmeV are well within the range 
where the resonance does exist, i.e., the smallness of the 
total spectral weight of the resonance does not matter. 

The subject of the present communication is the anal- 
ysis of the possible role of the spin resonance in the 
observed changes between normal and superconducting 



state in the optical data at high frequencies, lu ^ 0.5 — 
leV. Recently, Molegraaf et al. reported the results 
of their ellipsometry measurements on optimally doped 
and underdoped Bi2Sr2CaCu20s+5 with Tc = 88K and 
Tc = 66K respectively |l| . They observed that in the nor- 
mal state, the in-plane plasma frequency increases as 
with decreasing T, however below Tc it increases faster 
such that the actual value of ujpi{T — 0) is larger than 
the extrapolation from the normal state. The effect is 
very small: at optimal doping Sujpi ^ 10cto~^ is only 
1.3 X 10~^ of the plasma frequency ujpi ^ 7600cm~^ « 
leV, but detectable by the ellipsometry technique. 

The change of the plasma frequency in a superconduc- 
tor is related to the change of the fermionic self-energy 
between superconducting and normal states (see Eq. Q 
below). Conventional wisdom holds that at a; ujpi, 
which well exceed the magnetic bandwidth, interaction 
with low-energy spin fluctuations is not the dominant 
mechanism for the fermionic self-energy. We argue, how- 
ever, that while this is generally true for S(w) itself, 
Sujpi scales with S'S{uj) = Ssc(^^) — S„(ct;), and the latter 
comes from frequencies comparable to the superconduct- 
ing gap and can be captured within the low-energy, spin- 
fluctuation theory. We will see, however, that at high lu, 
the concept that intermediate fermions and bosons have 
equal energies fails, and a fermion with energy leV inter- 
acts with the whole band of magnetic fluctuations. As a 
result, STi^uj) scales with the integrated magnetic spectral 
weight, and the smallness of the spectral weight transfer 
into the resonance becomes crucial. We argue that the 
value of the observed shift of the plasma frequency is 
consistent with the fact that only 1 — 2% of the magnetic 
spectral weight is transferred into the resonance. 

Our reasoning is the following. At plasma frequency, 
the real part of the dielectric function e(w) changes sign. 
The dielectric function obeys e{Lu) = e(oo) -I- 47rz(T(ct;)/a;, 
where (t{lu) is the optical conductivity. By Kubo formula, 
a{uj) = ((w^,)V(47r))i?e[n(w)/(-iw)] where n(w) is fully 
renormalized current-current polarization operator and 
(Wpj)^ = Airne'^/m is the bare plasma frequency. The 
plasma frequency is then the solution of 



copi = jReIl{upi) (1) 
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To zero-order approximation, Il{ujpi) = 1, i.e., iVpi = 
Lu^i/ \/ e(oo). However, at any finite frequency \ — ReX\{uj) 
is still finite, and hence lo^i is sensitive to the change of 
the polarization operator upon entering the supercon- 
ducting state. This change of lo^i is small as supercon- 
ductivity mostly affects the form of n(a;) at frequencies 
comparable to the superconducting gap A ^ 0.04wpi. 

At high frequencies, w ~ leT^, normal and anomalous 
fermionic self-energies I](k, lj) « and $(k, cj) are 

both small compared to and to the leading order in 
the self-energy, the current-current correlator both in the 
normal and superconducting states is given by 

Jo S2 -I- S(ci^) -f - w) 

Substituting this expression into we find that 

the superconductivity-induced change in the plasma 
frequency at T = is related to the difference 
between superconducting and normal fermionic self- 
energies = - as 

^i^pi _ 1 
LOpl 2n^(cjp;) 

xRe / duj '- 2 (3) 

Jo [ujpi + ujZ^ + {u;pi - uj)Z^^^^^] 

where Sujpi — iv^'f — uj^i, li'ni^) is the real part of the 
polarization operator in the normal state, and = 1 + 
I]„(aj)/w is inverse quasiparticle residue in the normal 
state. By all accounts, at w copi, T,{lu) ^ uj, i.e., 
Z^ ^ 1. Hence, once the integral in © is dominated 
by frequencies where either lo or ujpi — lo are near ujpi 
(as we later verify), the precise form of the fermionic 
Z{ll!) at high frequencies does not matter, and StOpi/ujpi « 
SY^{ujpi)/(2ujpi). 

The computation of 6ujpi therefore reduces to the com- 
putation of the self-energy difference between normal and 
superconducting states. We present the result for ST,{uj) 
now and discuss its derivation later. We found that, 
within the spin-fluctuation scenario, there are two dis- 
tinct frequency regimes depending, roughly, on whether 
or not UJ exceeds the magnetic bandwidth. At small 
frequencies, the Eliashberg approximation is valid, and 
ST,{uj) is positive (see Fig. O. In this regime, internal 
fermions and bosons in the self-energy diagram have com- 
parable energies, i.e., a fermion is interacting only with 
spin fluctuations very near Q. We verified that if this be- 
havior extended up to a; ^ LUpi, Sujpi would be negative, 
in disagreement with the data. 

At larger frequencies, however, Eliashberg approxima- 
tion becomes invalid, and a novel (anti-Eliashberg) ap- 
proximation has to be used. In this approximation, we 
obtained that the change of the fermionic self-energy is 
proportional to the integrated change of the dynamical 
spin structure factor 

" / ^ 5: (4) 
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where g is the spin-fermion coupling constant estimated 
to be 5 - 0.7eV |l|l[ll|, and 5S,{c[,n) = 5f (q,ri) - 
^"(q, il) is the change in the dynamical structure factor 
5(q, n) ^ x"(q, + coth(w/2T)) in even {i = 1) and 
odd {i = 2) channels. The factor F{q) decreases away 
from Q but can be safely approximated by F{q) = 1 in 
the narrow momentum range where the resonance is ex- 
perimentally detectable. The integrated magnetic spec- 
tral weight near (tt, tt) is larger in the superconducting 
state 0,13, hence at high frequencies S'(w) is smaller in 
a superconductor. Using this JS(aj) we find that Sujpi is 
positive, in agreement with Q]. 

The crossover between the two regimes occurs at opti- 
mal doping at UJ = ujq ^ SOOmeT^, see Fig. |2| Theoreti- 
cally, we use one-band model to get Eq. ^ i.e., we assume 
that there is a frequency range above wq, where anti- 
Eliashberg approximation is valid, but interband transi- 
tions still can be neglected. The applicability of this ap- 
proximation has to be verified by comparing the results 
of one-band analysis with the data. 

The momentum and frequency integral in the r.h.s. 
of Q yields BS{S + l)/3 where B is the percentage of 
the spectral weight redistributed below Tc. Only the odd 
channel contributes to 5S{q, ft) and 1 — 2% of the spectral 
weight from this channel is redistributed p|, i.e., B ~ 
0.005 - 0.01. Substituting (5E(w) « -Sg^B/iiui) into 
the expression for ScOpi we obtain 5ujpi/ujpi ~ (1 — 2) x 
10"'^, in near perfect agreement with the experimental 
result 1.3 X 10"'^ P]. We emphasize that the agreement 
is entirely due to the fact that the integrated weight of 
the resonance is very small. If it wasn't, the blue shift 
of the plasma frequency would be much larger. To verify 
the prefactor, we went beyond estimates and evaluated 
the full integral in ^ using the normal state expression 
for Z^^: Z^ Kl + {iui/ujf^ where w ~ 0.35g. We 
obtained almost the same result for diOpi as above - the 
extra prefactor is 0.97 « 1. 

In Fig^ we plot the temperature dependence of 
5LJpi{T) together with the temperature dependence of the 
resonance peak intensity Q . The data for the fully inte- 
grated intensity are available for fewer temperatures and 
show roughly the same T dependence . We see that the 
T dependences of 5ujpi and of the resonance peak follow 
each other as it should be according to Eq. Q . 

We also searched for the explanation of the temper- 
ature dependences of Wp/ found by Molegraaf et al. in the 
normal state. The analysis of the neutron data shows 
that the dominant source of the temperature dependence 
of the spectral function is the T dependence of x"('^), 
while the one from (1 -I- cothti;/(2T)) is much weaker. 
Indeed, the neutron data presented for uo — 65meV 
show that the local x"(^) decreases by more than the fac- 
tor of 2 between 100 and 200K, when the thermal factor 
coth(a;/2T) w 1 for both temperatures. It is therefore 
likely that the T dependence in the normal state orig- 
inates in thermal corrections to the parameters of the 
dynamical spin susceptibility (i.e., to the magnetic cor- 
relation length). Theoretically, these corrections form 
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FIG. 1: The temperature dependence of the 

superconductivity-induced change of the plasma frequency 
in _Bi2212 M vs the change in the resonance peak intensity 
for YBCO Upper panel - optimal doping, lower panel - 
underdoped samples. We plot YBCO data for two different 
underdoped samples for comparison. 



regular series in {T/EpY [1|> hence the temperature de- 
pendence of ojpi must be T^, unless superconductivity 
interferes. 

We now turn to the conductivity. Molegraaf et al. re- 
ported that the Drude optical spectral weight integrated 
up to a frequency of lOOOOcm"^ (including the conden- 
sate contribution) increases below Tc, and argued that 
this increase is compensated by the decrease of the op- 
tical spectral weight integrated between lOOOOcm"^ and 
20000cTO~^, where interband transitions are overly rel- 
evant. This, they argued, implies that the condensate 
contribution to cri(w) is not compensated by the reduc- 
tion of the conductivity in the superconducting state up 
to frequencies when one-band description fails. 

Our results are inconsistent with this claim. The opti- 
cal integral 5I{Vt) = h? Sa{Lu)du! for ~ leV covers 
the range of frequencies where the Eliashberg approxima- 
tion is valid (SIi) , and the range where Eq. Q is valid 
((5/2). Using Eq. we can compute the partial inte- 
gral 6l2{^) = Tt? 6u(ijj)dijj where, we remind, loq is the 
lower boundary for the applicability of Eq. . The ap- 
proximations for a are less robust than for ujpi as for the 
conductivity we need to know ImZ^ at high frequencies, 
where the use of the low-energy theory is questionable. 
Still, at high frequencies, Im[l/ Z{uo)] k, —ImT,{uj)/uj is 
negative, hence, according to the high-frequency con- 
ductivity a{Lo) ^ 1/lu'^t{lu) oc is larger in the su- 
perconducting state than in the normal state extrapo- 
lated to T = 0. Accordingly, dl2{^) is positive. Sub- 
stituting Eqs. l|2Jl and 10} into the Kubo formula, and 
evaluating the integrals with the same Z(io) as above, we 
obtain Shifl) « (2 - 4) x IQ-^ieVf for n ~ leV. 

The integral for (5/2 (fi) converges at the upper limit; 
formally taking U — 00 yields almost the same dl2- As 
the /—sum rule must be satisfied within the one-band 
model (if the bandwidth is formally set to infinity), the 
value —612(00) gives the estimate of the optical inte- 
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FIG. 2: Right panel - the difference in the electronic self- 
energy between normal and superconducting states calculated 
explicitly and using Eliashberg and anti-Eliashberg approxi- 
mations. The frequency is measured in units of lD ~ 0.35g ~ 
250meV. Clearly, the anti-Eliashberg approximation is much 
better at large frequencies. Left panel - the low-frequency 
region, where the Eliashberg approximation is valid. 



gral SIi{u!o) over frequencies where the Eliashberg theory 
is valid. From our consideration, this integral is nega- 
tive, i.e., the condensate contribution is overcompensated 
by the reduction of the conductivity in the supercon- 
ducting state below ujq 0|. As an independent check, 
we computed optical integral 6I1 within Eliashberg the- 
ory and indeed obtained that SIi changes sign at 
uj ~ 300meV at optimal doping (see FiglJ)). The ac- 
curacy of our Eliashberg calculations is not sufficient to 
compare the values of the two contributions to 61. 

Our results therefore indicate that at T = 0, the full 
optical integral is exchausted at frequencies smaller than 
the bandwidth. The optical integral evaluated over larger 
frequencies is larger at T = than in the normal state 
extrapolated to T = 0. The crossover frequency luq some- 
what increases with underdoping (see caption for Fig|31) 
but theoretically it still remains 0((D) even in strongly 
underdoped materials Note that this does not con- 
tradict the idea that that superconductivity is driven by 
the decrease of the kinetic energy ^Gj , as within the same 
approach, the decrease of the kinetic energy also comes 
from frequencies (1 — 2)w [l7| . 

We now describe the calculation of the self-energy, Eq. 
Q in some detail. We assume that the fermionic self- 
energy predominantly comes from the fermion-fermion 
interaction in the spin channel and can be viewed as 
being mediated by spin collective modes with momenta 
near (tt, tt). Quite generally, the imaginary part of the 
fermionic self-energy is given by 

S"(k,c.) = ^ Jg\"{q,n)G"{k + ci,u: + n) 

dQ(fq (5) 



tanh ■ 



n 



2T 



coth — 
2T 



where g is the fully renormalized vertex, xCq? ^) is the 
propagator of the collective mode, and G(k + q,uj + il) 
is the full fermionic Green's function. 

At small frequencies, the physics that allows one to ne- 
glect vertex corrections (i.e., set g = g), is the fact that 
at strong coupling, spin fluctuations are overdamped in 
the normal state and are slow modes compared to elec- 
trons. Hence, an effective Migdal theorem is valid |^. 
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FIG. 3: The results of Eliashberg calculations for coupling 
A = 1. Left panel - self-energy, right panel - the optical inte- 
gral 51{uj) (including the condensate piece). For A — 2, SI{ui) 
changes sign at a larger ui ~ 2uj. 



By the same reason, the momentum integration in Eqn 
(O is factorized - the integration transverse to the Fermi 
surface involves only fast fermion, while the integration 
along the Fermi surface is over a slow bosonic momen- 
tum. This computational procedure is called Eliashberg 
approximation. It describes the jump in S" at A -I- ujres 
that is the key element of the peak/dip/hump behav- 
ior 0,1^. However, for S'E,{uj) this approximation is valid 
only as long as external fermionic frequency u is smaller 
than a typical frequency at which the momentum inte- 
gral of Sx"{^) converges. At strong coupling, this gen- 
erally happens at Q comparable to the effective bosonic 
bandwidth Q defined such that in the normal state, 
becomes less than uj at lo > u), i,e., spin-fluctuation scat- 
tering becomes ineffective. When spin-fermion coupling 
g is less than the fermionic bandwidth VF, (D ~ 0.35g, in 
the opposite limit a) ~ W^/g ~ J. Note that within the 
same model, A ~ O.Sw Q. 

When a; > w, and k « kp, G"(k + q,uj + n) w 
Im{ujZ^)~^ , and it can be taken out of the integral. To 
the same accuracy, the T dependent factor in (|SJ| is ap- 
proximated by 1 + coth(ri/2r). We then obtain: 



Im [(5S(w)] 
dQd^q 



-3g Im 



87r3 



(6) 



where F(q) subject to i^(Q) = 1 decreases at |Q — q| = 
O(IQI) and reflects the fact that the spin-fermion model 



is only valid for bosonic momenta near (tt, tt). As the 
Kramers- Kronig transform of © is infrared convergent, 
the full self-energy is given by 

We find therefore that at high frequencies w ^ tD, the 
correct computational procedure for S'E,{u;) is opposite 
to the Eliashberg approximation ~ instead of factorizing 
the momentum integral, one can neglect the momentum 
dependence in the Green's function and perform the full 
2D momentum integration over the bosonic momenta. 
We verified that in this, anti-Eliashberg approximation, 
vertex corrections are again small, this time in lo/lo, such 
that g = g in Eq. © . 

Obviously, there should be a crossover between Eliash- 
berg and anti-Eliashberg approximations as frequency 
increases. To understand where it is located, we eval- 
uated ST,"{uj) explicitly, using the normal and supercon- 
ducting forms of the dynamical spin susceptibility ob- 
tained earlier and compared the full result with the 
two approximate forms. The results are presented in 
Fig. 13 We see that for ui < uj, Eliashberg approxi- 
mation is much closer to the full result. However, for 
uj > {1 — 2)u) « 250 — 500meV the Eliashberg approxima- 
tion is well off, while the anti-Eliashberg approximation 
is rather close to the full expression. This justifies our 
use of Eqn. Q for optical properties above 500meV. 

To summarize, in this paper we argued that the 
superconductivity-induced blue shift of the plasma fre- 
quency, detected in the ellipsometry studies, can be ex- 
plained within the magnetic scenario for the cuprates. 
We found that Sujpi scales with the change of the in- 
tegrated magnetic spectral weight (5S'(q, fi). The mag- 
nitude of Sujpi/ojpi is small {r^ 10"'^) as the integrated 
SS{q, ft) accounts for only a small fraction of the total 
spectral weight. We also predict that the optical integral 
converges below 0.5eV. Careful measurements of this in- 
tegral should either confirm or disprove our claim. 
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